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Soft-photon corrections in multi-body meson decays 
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\Q \ Abstract 

The effects due to soft-photon emission (and the related virtual corrections) in multi- 
body decays of B, D, and K mesons are analysed. We present analytic expressions 
for the universal 0(a) correction factors which can be applied to all multi-body decay 
modes where a tight soft-photon energy cut in the decaying-particle rest-frame is applied. 
All-order resummations valid in the limit of small and large velocities of the final-state 
particles are also discussed. The phenomenological implications of these correction factors 
in the distortion of Dalitz-plot distributions of K — > 37r decays are briefly analysed. 
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In the last few years the large amount of data collected at flavour factories has allowed to 
reach statistical accuracies around or below the percent level in several decays modes of B, 
D, and K mesons. At this level of accuracy electromagnetic effects cannot be neglected. 
The theoretical evaluation of these effects is a key ingredient to extract from data a precise 
information about weak interactions or strong dynamics, such as the determination of 
CKM matrix elements or the extraction of of tttt scattering lengths. 

The theoretical treatment of the infrared singularities generated within QED is a well 
known subject and one of the pillars of quantum field theory. A clear and very general 
discussion can be found, for instance, in Ref. [1,2]. These general properties of QED 
have been exploited in great detail in the case of genuine electroweak processes, or pro- 
cesses which can be fully described within perturbation theory within the Standard Model 
(SM). More recently, a similar program has been extended to a few decay modes of K 
and B mesons (see e.g. Ref. [3-7]), which can be described within appropriate effective 
field theories (EFT). The purpose of the present article is to complement and gener- 
alise these EFT studies, analysing the general structure of electromagnetic corrections in 
multi-particle final states. In particular, we are interested in the distortions of the non- 
radiative decay distributions (Dalitz plot parameters, form factor slopes, etc. . . ) induced 
by electromagnetic effects. To a large extent, these effects have a universal (long-distance) 
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character: their structure can be evaluated independently of the short- distance dynamics 
which originate the meson decay 



2 The photon-inclusive decay distribution at 0(a) 

From the experimental point of view, the most convenient infrared-safe observable related 
to the process P — > Pi . . . P/v is the differential photon-inclusive distribution 

dr incl ( SlJ ; P max ) = dr(P - P x . . . P N + n 7 ) | E £ 7<£m a X , (1) 

namely the differential width for the process Po — > P\ . . . P/v accompanied by any number 
of (undetected) photons, with total missing energy less or equal to E max in the Po rest 
frame. In addition to E ma,x , the differential photon-inclusive distribution depends on 
kinematical variables describing the visible particles. A convenient choice for the latter 
ia3 

ij l(Po-Pi) 2 i = 0,^0. [) 

The photon-inclusive distribution in can be decomposed as the product of two 
theoretical quantities: the so-called non-radiative width, dr°(sjj), which survives in the 
a — > limit, and the corresponding energy-dependent electromagnetic correction factor 
E 

At any order in the perturbative expansion in a the energy dependence of Q(sij] E) is 
unambiguous and universal up to terms which vanish in the limit E — > [2]. The P- 
independent part of Q contains both universal terms, such as the Coulomb corrections, 
and non-universal terms depending on the short-distance dynamics which originate the 
decay. In order to discuss the separation between universal and non-universal terms, 
we start presenting the calculation of Q(sij]E) at 0(a) in the limit of a real point-like 
effective weak vertex. 

The general decomposition of Q(sij\ E) at 0(a) is 



dr incl (^; p max ) = dr°( Sii ) x n( Sij] p max ) . (3) 
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Q{ Sij ;E) = l + J2QiQiM s if E ), ( 4 ) 



i,j=0 



M** ; = I [ 2 % ln (S) + F « + Hg + H<} + H™ + (E) ] , (5) 



where Qi^o are the charges of the final-state particles in units of e and Qo = — J2iLi Qi- 
The terms 6y and Fij are unambiguously determined by the real-emission amplitude, 
while the Py functions are associated to virtual corrections. 



1 Although redundant, this choice of variables allows us to keep the discussion on a general ground. 
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The S matrix element corresponding to the emission of a real photon can be decom- 
posed as 



N 



out (Pi(pi) • • • P N (Pn) + 7(e, k)\P (p )) in = -ieM xKx (27r) 4 <5 4 (p - £> ~ k ) ' ( 6 ) 

i=i 

where .Mo is the invariant amplitude of the non-radiative process: 

N 

out (P 1 ( Pl ) . . . Piv(Pn)|^o(Po))in = -^0 (2vr) 4 5 4 ( Po " J>i) 

and 



i=i 



(7) 
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(8) 



The integration of the real-emission amplitude in the soft-photon approximation with 
a photon-energy cut E (namely neglecting O(E) terms) and regularizing the infrared- 
singularities with a photon mass m 7 , leads to 



dT rca \s l:j ;E) = dT°(s t3 )x [ 

Je 



d 3 k 



, <E (2vr)3 2E, 



£1*1 



spins 
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dr °(%) >< £ E K ln (il) + ^ + <w 



(9) 



ij=0 



where [1,2] 
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4m 2 m 2 



( Sij - m? - m 2 )' 



1/2 
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The finite term Fjj depends on the specific cut applied on the (soft) photon energy. 
Imposing the condition p Q ■ k < m E, corresponding to a cut in the Pq rest frame, leads 
to 



A,,- / dz 



e(z) 



In 



x p(z)[e 2 (^) -p 2 (z)] V e (^)-p(^)y 



e(z) +p(z) 



where 

e(*) = 
p(z) = 
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[7w(l - z) +7oj(l + z)] , 



;i-/%) 1/2 ' 



(ii) 
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K-l)(l-^) 2 + (7o 2 ,-l)(l + ^) 2 ] 

1/2 



+2 ( 70,70,^^ - A 4j ) (1 + z)(l - z) 



"i j 



— m} — m 2 
2s 



(13) 



with the special case i = j given by 
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As far as virtual corrections are concerned, the universal infrared singular term cancels 
out the ln(m 7 ) dependence in Eq. ([9]), and the remaining finite terms are encoded into 
the three functions in Eq. (jSJ). Regularizing UV divergences by means of dimensional 
regularization and renormalizing the real point-like weak vertex in the MS scheme leads 
to 
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where 
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ln(l - £) 
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H^, which is singular in the limit of vanishing velocity among any 



The first term 

pair of charged particles, is a genuine long-distance effect: it corresponds to the Coulomb 
interaction among the two charged particles. This term can indeed be evaluated also in 
non-relativistic quantum mechanics by means of semi-classical methods (see e.g. Ref. [8]). 

The second term, H^j , which depends explicitly on the ultraviolet renormalization 
scale /i, is manifestly not universal: its scale dependence cancels out in Eq. ([3]), or in the 
physical observable, by the corresponding scale dependence of the weak amplitude. The 
finite 0{ot) term resulting after this cancellation cannot be computed without knowing 
the short-distance behaviour of the amplitude. Note that, in the approximation of a 
point-like weak vertex, this missing piece affects only the overall normalization of the 
photon-inclusive distribution and not its kinematical structure. 

By construction, H}^ is what remains after isolating the manifestly universal and 
manifestly non-universal terms H$ and . More explicitly, Hj^ is the finite part 
of the universal three-point function function after subtracting ultraviolet and infrared 
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divergences and the Coulomb termH 

d d k 



H, 



IK 




(2 Pi + k)^2 Pj - ky 



i(2n) d [(pi + k) 2 - mf] [(pj - k) 2 - m 2 ] [k 2 - m 2 ] 



m 



(18) 



3 Resummations and universal correction factor 

The E — > singular terms in Eq. (jSJ) and the ftij — *> singular terms in H£, which 
represent the potentially largest correction factors, can be summed to all orders in a. 

As shown in Ref. [2], the resummation of the a n \n n (E) terms allows us to remove the 
E — > singularity, giving rise to the following exponential term 



(2E 

Q B (^;£) = — 
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i,j=0 



(19) 



The resummation of the (na/j3ij) n Coulomb terms is encoded by the semi-classical 
result [81 
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The Pij — > singularity does not disappear and it is strengthened in the case of opposite- 
sign charges (attractive interaction), but it remains an integrable singularity over the 
final-state phase space. 

The two resummed expressions in Eqs. (fT9j) and (120]) are relevant in two different 
kinematical regimes: Qc(%) is relevant in the /3y — > limit, while QB{sij', E) acquires a 
non-trivial kinematical dependence only in the — > 1 limit. We can therefore factor- 
ize the two effects up to sub-leading 0(a 2 ) corrections. This allows us to consider the 
following generalization for the universal part of the electromagnetic correction factor: 



Q eS (sij;E) = Q B (sij;E) x Oc(sij) x 
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l + -Y,QiQj{Fij + H, 

71 * — ' 

i,3=0 



IR\ 

ij ) 



(21) 



This expression provides a good description of the leading kinematical corrections 
induced by soft photons in multi-body meson decay. The approximations/validity-limits 
of fl e s(sij] E) can be listed as follows: 

• The leading kinematical singularities, namely the a n / terms for — > and the 

1, are summed to all orders. 



a n ln n (E/m ) n ln n (l - /%) terms for 
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2 The result in Eq. (|17|l is is valid only for Sy variables in the physical range, namely s%j real and 
positive (such that all terms in Eq. (fl7|) are real). 
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• The regular contribution of the real photon emission (Fy) is correct up to con- 
stant terms of 0(a 2 ) and energy- dependent terms of 0(aE/A), where A is a typical 
hadronic scale. More precisely, the corrections linear in E are controlled by the 
derivatives of the non-radiative amplitude with respect the kinematical variables: 
0(aE x dA/dsi) [9]. In several cases the tightness on the photon-energy cut nec- 
essary to keep these corrections under control can thus be quantitatively controlled 
by the smoothness of the non-radiative amplitude. In practice, the photon-energy 
cut is rarely a problem in 7r and K decaysj, while it is a non-trivial constraint for 
heavier mesons. 

• The virtual corrections encoded in H\J- are only the universal contribution of low- 
energy photons within an effective theory valid below the scale A (A < M p ), with 
real effective couplings in the a — > limit. High-energy modes provides in gen- 
eral additional infrared-safe 0(a) contributions which should be evaluated mode 
by mode (non-universal terms), and which are different in case of final-state lep- 
tons or mesons@ By an appropriate matching procedure, these additional terms 
can be reabsorbed into the normalization and the kinematical dependence of the 
non-radiative amplitude. In light meson (n and K) decays these extra terms are 
necessarily smooth functions of the kinematical variables 0(aSi/A 2 ) and thus can 
be safely neglected. These ultraviolet effects are potentially larger in heavy me- 
son decays, but also in this case they are subleading with respect to the leading 
logarithmic singularities included in fl e g(sij; E). 

• The only cases where virtual effects not included in Eq. (j2"T|) are potentially relevant 
are the singular points corresponding to the formation of Coulomb bound states. A 
notable example is the pionium formation [10], which has recently been observed 
in K — > 3tc decays [11,12]. Such states are treated here as different final states, 
which should be eliminated by appropriate kinematical cuts (as done for instance 
in Ref. [11,12]). Given the extremely narrow widths of Coulomb bound states, and 
the low probability formation, these effects are relevant only in very tiny regions of 
the space space and can be safely neglected in heavy-meson decays. 

4 A specific application: K + —> tt + 7t + 7t~ decays 

The high-statistics measurements of the K — > 37r Dalitz Plot distributions performed by 
the NA48/2 collaboration [11,12] have recently received a considerable attention because 
of the possibility to extract a precise information on tttt scattering lengths [13-17]. 

3 The only exceptions are modes where the bremsstrahlung is strongly suppressed compared to the the 
direct emission by symmetry arguments, such as the helicity-suppressed K —f ev^) or the CP-violating 

Kl — > 7T + 7T~(7). 

4 Having assumed real effective couplings in the a — > limit, we have also ignored the O(a) electro- 
magnetic corrections to the strong phases of the amplitude. For smooth strong phases these can be easily 
be incorporated starting from the imaginary part of the three-point function in Eq. (|18p , as discussed for 
instance in Ref. [3] for the K — > tttt case. 
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Figure 1: Radiative corrections in the K 




7r + 7r + 7r decay. Left: density plot of 



[fi e ff(%; E) — 1], evaluated with the full correction term in Eq. (121]) with E = 5 MeV. 
Right: density plot of [Qc(sij) — 1] (Coulomb term only). The s« are in units of uik- 



The leading mechanism which allow to measure irn scattering lengths (and particularly 
the ao — a,2 combination) in K — > 37r decays is the 7r + 7r~ — > 7r°7r° re-scattering at the 
7T + 7r _ threshold, which produces a prominent cusp in the M n o n o distribution of the K + — > 
7i + 7i°7i° decay [13]. The strength of this singularity is proportional to ao — 02, but also 
to phenomenological parameters introduced to describe the K + — > 7r + 7r + 7r~ amplitude 
(see Ref. [14, 15]). The latter must be determined by experiments from a fit to the 
K + — > 7r + 7r + 7r _ decay distribution, which is likely to receive sizable electromagnetic 
distortions because of the three charged particles in the final state. 

In Fig. [TJwe show the impact of soft-photon corrections in the K~ — > 7r + 7r + 7r~ decay 
distribution. In particular, we compare the result obtained with the full universal cor- 
rections factor in Eq. (T2T]) or using only the Coulomb term in Eq. (1201) . As expected, 
radiative corrections induce sizable distortions, especially at the border of the Dalitz 
plot distribution. However, these are well described by the Coulomb term up an overall 
normalization factor of (9(1%). The procedure adopted by the NA8/2 Collaboration to 
correct K + — > 7r + 7r + 7r~ data using only the Coulomb term is therefore well justified a 
posteriori. 

As discussed in the previous section, our general treatment do note take into ac- 
count the formation of Coulomb bound states. Such process occur at the border of the 
K + — > 7r + 7r + 7r~ Dalitz plot, when one of the two 7r + 7r~ pairs is at rest. In order to 
determine the K + — > 7t + 7t + tt~ decay parameters relevant to the analysis of Ref. [14,15], 
the narrow regions at the border of the Dalitz plot with Coulomb corrections of (9(100%) 
should therefore be eliminated by appropriate kinematical cuts. This procedure is per- 
fectly consistent with the cut of the pionium region (around the peak of the M n o n o cusp) 
performed in Ref. [11]. 
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